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AERONAUTIC  SYMBOLS 
1.  FUNDAMENTAL  and  derived  units 


.•1 


Symbol 

Metric 

English  y  . 

Unit 

Abbrevia¬ 

tion 

Unit 

Abbrevia- 

tion 

Length.. ...A 

Time.... . 

Force _ 

;  "ii  ;  ■ 

t  .. 

F  . 

meter.  --  -----  - 

second - — - 

weight  of  1  kilogram..,— 

m 

s 

hg  . 

•'  ■- 

foot  (or  mile) _ - _ 

second  (or  hour) . . 

,  weight  of  1  pound _ 

ft,  (or^f;)V 
sec.  (or  hr.) 

ib.," 

-  yy 

h  Power^ 
Speed.— Li— 

.  \  “  \  v  ■ 

V*  ; 

■  .  .  K 

horsepower  (metric)  -—- 
/kilometers  per  hour.— -- 
\meters  per  second - — 

k.p.h. 

m.p.s. 

horsepower _ 

miles  per  hour.. _ I 

feet  per  second . . 

lO:‘-  ■  1 

**r‘ 

■  n>.p.b^  .  y 

•  ^  t.'-.v -r 

2i  GENERAL  SYMBOLS 

v,  Kinematic  viscosity 


i 


g  Specific  weight  of  “standard  ”  air,  1.2255  kg/m*  or 

'  is  of  : %f:  0.0765llb.Aaiit.»,  ,  *  : v 


Mi 

'si 

S>w> 

G, 

■b, 

c, 

V 

£' 

V, 

2> 

L, 

I>, 

A, 

Dp, 

C, 

R, 


•Moment  of  |n^a,-,inP.  j(Indicaie  oris  oi 
'  radius  of  gyration  k  by  proper  subscript.) 

:  ,S  Coefficient  of  viscosity 


•94': 


:•  v\: 


Area 

Area  of  wing 

Gap 

v  Span 1  ■/SV'Vv- U  .vy?:':.-  ’  "V?- :  'U  'l  \ 

Chord 

Aspect  ratio  , 

True  air  speed  r.'-y-y ,y  y' 

Dynamic  pressure  —  ypV* 

Lift,  absolute  coefficient 

'  s'  .  •f'  y  ■  <n  ■  v.  •• 

Drag,  absolute  coefficient  CD  — 

Profile  drag,  absolute  coefficient  —  ^ 
Induced  drag,  absolute  coefficient 

^T>  />  ^ 


t:1-?  ;v 

-i4i>  aerodynamic  symbols 


i„,  Angle  of  setting  of  ’.rings  (relative  to  thrust 

V  line)  •  '  ;  ,  ft 


Q, 

*V 

vi 

V 


Uv 'Jo- 

.. ;  -v.; ' 


Parasite  drag,  absolute  coefficient  CD  — 


Cross-wind  force,  absolute  coefficient  Gcx 
Resultant  force 


C 


<*> 

Of0; 

Off, 

an} 

7, 


Angle  of  stabilizer  setting  (relative  to  thrust 
line) 

Resultant  moment 
Resultant  angular  velocity  f '  ; 

Reynolds  Number,  where  l  is  a  linear  dimension 
(e.g.,  for  a  model  airfoil  3  lif.t^horii,  100 
m.p.h,  normal  pressure  at.  15°  .Cf,  the  cor- 
1  responding  number  is  234,00b;  ^  j6r  a  model 
of  10  cm  chord,  40  m.p.s.  the  corresponding 
..  number  is  274,000)  ‘  x 

Center-of-pressurc  coefficient  (ratio  of  distance 
of  cvp.from  leading  edge  to  chord  length) 
Angle  of  attack  ‘  ^ 

Angie  of  down  wash  v 

Angle  of  attack,  infinite  aspect  ratib  v  '  r  •;  / 
'Angle  of  attack,  induced  ’  -AH- 

Angle  of  attack,  absolute  (measured  from  zero- 
lift  position)  .  .  v..;v 

Flight-path  angle  ^  ',.7  v 
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AN  ANALYSIS  OF  LONGITUDINAL  STABILITY  IN  POWER-OFF  FLIGHT 
WITH  CHARTS  FOR  USE  IN  DESIGN 

By  Charles  H.  Zimmerman 


SUMMARY 

This  report  presents  a  discussion  of  longitudinal 
stability  in  gliding  flight  together  with  a  series  of  charts 
with  which  the  stability  characteristics  of  any  airplane 
may  be  readily  estimated. 

The  first  portion  of  this  report  is  intended ,  primarily , 
for  students  of  the  subject.  The  relationships  governing 
stability  characteristics  are  derived  from  equations  of 
equilibrium  referred  to  moving  axes  that  are  tangent 
and  perpendicular  to  the  instantaneous  flight  path.  It 
is  shown  that  instability  of  the  motion  can  arise  only 
through  an  increase  of  linear  and  angular  momentum  in 
the  system  during  one  complete  cycle.  The  interaction 
of  events  leading  to  increase  or  decrease  of  momentum 
during  a  cycle  is  explained  in  detail.  The  construction 
of  charts  showing  the  effects  of  the  nondimensional 

parameters  CL,  CD,  —  mq,  and  —\xma  upon 

the  stability  characteristics  is  explained  and  the  effects 
of  the  more  important  of  the  aerodynamic  and  mass 
characteristics  of  the  airplane,  as  revealed  by  the  charts, 
are  discussed. 

The  latter  portion  of  the  paper  is  devoted  to  a  series  of 
40  related  charts  with  which  the  dynamic  stability  of  any 
airplane  in  power-off  flight  may  be  readily  estimated. 
The  use  of  the  charts  is  explained  in  detail  so  that 
reference  to  the  earlier  discussion  is  unnecessary. 

INTRODUCTION 

The  longitudinal  stability  of  aircraft  has  received 
very  extensive  and  exhaustive  treatment  by  able 
writers  (see  references  and  bibliography),  but  the 
classical  treatment  of  the  subject  has  been  rather 
difficult  for  those  not  familiar  witli  higher  mathematics. 
The  study  reported  herein  was  undertaken  with  the 
purpose  of  making  more  understandable  the  mathe¬ 
matical  treatment  and  of  preparing  a  method  of  esti¬ 
mating  stability  characteristics  that  would  be  suffi¬ 
ciently  accurate  and  rapid  to  appeal  to  practical 
designers. 

The  section  preceding  the  group  of  charts  for  deter¬ 
mining  stability  characteristics  in  power-off  flight  is 
devoted  to  a  derivation,  in  relatively  simple  terms, 
of  the  mathematical  relationships  and  to  a  discussion 
of  the  formulas.  The  portion  following  the  group  of 
charts  consists  of  an  explanation  of  the  method  of 
using  them.  It  is  not  necessary  to  read  the  first  por¬ 
tion  in  order  to  use  the  charts  with  satisfactory  results. 


All  symbols  not  given  in  the  report  cover  are  defined 
where  used  and  are  also  listed  in  the  appendix. 

I.  ANALYSIS  AND  DISCUSSION 

DERIVATION  OF  MATHEMATICAL  FORMULAS 

Definition  of  stability  characteristics, — The  sta¬ 
bility  characteristics  of  an  airplane  are  those  qualities 
which  define  the  nature  of  the  motion  after  a  deviation 
from  an  initial  condition  of  equilibrium.  The  motion 
may  be  periodic,  consisting  of  a  series  of  oscillations 
having  a  certain  period  and  rate  of  increase  or  decrease 
in  amplitude,  or  aperiodic  with  a  certain  rate  of  return 
toward  or  deviation  from  the  equilibrium  position. 
In  many  stable  airplanes  the  return  to  a  condition  of 
nonoscillating  equilibrium  is  spoken  of  as  aperiodic  or 
“dead  beat”  when  it  is  essentially  oscillatory  in 
character  but  very  heavily  damped. 

Fundamental  concepts  and  assumptions. — The  forces 
and  moments  determining  the  motion  of  the  airplane 
are  of  two  kinds:  (1)  Aerodynamic  forces  and  moments 
created  by  movement  of  the  lifting  and  control  sur¬ 
faces  relative  to  the  surrounding  air;  (2)  mass  forces 
and  moments  arising  from  the  weight  and  acceleration, 
angular  as  well  as  linear,  of  the  airplane.  The  funda¬ 
mental  basis  of  the  discussion  presented  in  this  report 
is  that  at  all  times  there  exists  a  state  of  equilibrium 
between  the  mass  forces  and  moments  and  the  aero¬ 
dynamic  forces  and  moments. 

A  complete  treatment  of  the  stability  of  airplanes 
would  be  extremely  lengthy  and  very  complex.  Cer¬ 
tain  assumptions  have  therefore  been  made.  As  the 
motion  of  an  airplane  is  three  dimensional,  it  is  to  be 
expected  that  any  treatment  of  the  subject  will  be 
incomplete  if  if  neglects  certain  of  the  components  of 
the  motion.  Fortunately,  conventional  airplanes  are 
symmetrical  (within  limits  here  applicable)  with 
respect  to  the  plane  that  includes  the  fuselage  axis  and 
is  perpendicular  to  the  span  axis.  It  is  obvious  that  a 
longitudinal  motion  having  no  component 'of  linear 
yelocity  perpendicular  to  that  plane  or  no  component 
of  angular  velocity  about  any  axis  lying  in  that  plane 
cannot  introduce  asymmetric  forces  or  moments. 
Such  motion  can  therefore  be  treated  as  an  independent 
phenomenon. 

The  longitudinal-stability  characteristics  will  neces¬ 
sarily  be  affected  by  any  deflection  of  the  lifting  or 
control  surfaces.  The  influence  of  wing  elasticity  and 
of  free  longitudinal  control  will  not  be  considered  in 
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the  primary  analysis  because  such  a  consideration 
would  complicate  the  relationships  and  obscure  the 
fundamental  principles. 

It  is  assumed  that  forces  and  moments  acting  upon 
the  wing  and  the  horizontal  tail  surfaces  vary  as  the 
square  of  the  air  speed  and  the  first  power  of  the  angle 
of  attack  of  the  individual  surfaces  and  that  they  are 
not  affected  by  the  rate  of  change  of  either  the  air 
speed  or  the  angle  of  attack.  Also  the  forces  upon 
the  lifting  surfaces  are  assumed  not  to  be  affected  by 
the  rate  of  rotation  of  those  surfaces  (reference  i). 

Each  of  the  foregoing  assumptions  necessarily  in¬ 
volves  a  certain  degree  of  approximation  but  they  are 
confirmed  by  comparison  between  measured  and  cal¬ 
culated  values  of  stability  characteristics  (reference  2 
and  unpublished  data)  and  are  justified  by  the  simpli¬ 
fication  of  the  relationships  they  permit. 

Equations  of  equilibrium. — As  has  been  previously 
stated,  the  course  of  the  airplane  in  flight  is  determined 
by  the  conditions  necessary  to  maintain  equilibrium 
between  mass  and  aerodynamic  forces  and  moments 
at  all  times.  In  steady  flight  the  equilibrium  may  be 
expressed  by  the  equations  (see  fig.  1): 


Figure  1.— Angular  and  vectorial  relationships  in  liight,  power  oil. 


W  sin  7+|p^C,„=0  n 

M'cos  7— b  (l) 
IpV’&O.-O  c 


where  (la)  refers  to  forces  tangent  to  the  instantane¬ 
ous  flight  path,  (lb)  to  forces  perpendicular  to  tlu 
instantaneous  flight  path  in  the  plane  of  symmetry 
and  (lc)  to  moments  about  an  axis  through  the  centei 
of  gravity  and  perpendicular  to  the  plane  of  symmetry 
After  displacement  from  the  steady  flight  condition 
the  equations  of  equilibrium  read: 

1  dV 

W  sin  (7  + AY^^P^d- Al^^CDd- ACd)^  m~dt  a 

\V cos (y Ay) — 7)p(V 4~  A V)2S(Ol 4" A C£)  —  mV ^ 

^p(V+AV)2Sc(C„+ACn)=mkr^p  Cj 

where 

acceleration  tangent  to  flight  path. 

ydy  centrifugal  acceleration  normal  to  the  flight 
V  dt'  path. 

d?B  angular  acceleration  of  airplane  about  tlu 
dt2’  lateral  axis. 

Since  the  effects  of  angular  velocity  and  of  acceleration 
upon  the  forces  are  neglected,  ACd  may  be  written  as 

Aa^and  A(7L  as  A Cm  may  be  written  as 

Ag^2  where  A 3=2,  the  angular  velocity  in 

pitch,  since  q  is  zero  in  the  original  condition;  and 
AV,  Ay,  Act,  and  q  are  small  quantities  by  assumption. 
Terms  involving  products  of  two  or  more  of  them  will 
be  neglected. 

Also 

sin  (7 -f  Ay)  =  sin  7  cos  Ay+cos  7  sin  A7 
=sin  7-j- A7  cos  7 

|  and 

cos  (7  + A7)  —  cos  7  cos  Ay  —  sin  7  sin  A7 
—  cos  7 — A7  sin  7 

I  Then: 


W  sin  y  +  WAy  cos  y  -\--^pV2SCD-r  pVSCDAV  -V^p\  2S-j^Aa—  m  ^ 

1  1  dCr  dy 

U’  cos  7—  IKA7  sin  7  —  9P  V  2SCl  —  p\  SClA\  —  2^  ^  da  A«—  “  m\  ^ 

\jV^ScCm+PVScCmAV+\pV^cd^Aa+\jViScd^q=mky-^ 


Subtracting  fl)  from  (3) 

XVAy  cos  y-\-pVSCDAV-\~~p\  ^-^Aa^—m-^  a 

-  WAy  sin  y-9VSCLAV-\f>V*d~±A*=  ~mVJt  b  (4) 

pVScCmAV+\p\^Scd^q+\prScd^Aa=mk/^  c 
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From  equations  (1) 


W  cos  y=-pV2SCL 
W  sin  7=  —-^pV^SCo 


sideration.  The  change  in  moment  due  to  the  change 
in  the  angle  of  attack  at  the  tail  is 


Therefore, 


-  ip V>SCL  Ay  -pVSCn  A  V—  IpPS— Aa = ml~ 
ipV*SCDAy-pVSCl.AV-ipV*S&£Aa=  -mVjjj  b  (5) 

Dividing  (5a)  and  (5b)  by  ^  pV^S  and  (oc)  by  pVS 


.ClAy-2CD^-d^Aa- 


dV 

9  m  (it 
pS?  V 


lMt=~p  VHS  tv  j^—Aa  t 

-  OCKj 

where  slope  of  normal-force  curve  of  tail  surfaces. 

dot  i 

Slf  area  of  horizontal  tail  surfaces. 

77 1,  tail  efficiency, 
or,  in  coefficient  form, 

A p  S,  dCzt  ql 

A Cm'~l  Sv‘d^7  V 

Introducing  an  empirical  factor  K  to  allow  for  wing 
damping 

dC„  Kl  .S',.  dCz,  l 
~d£-  c  SV‘da,  V 

Therefore  (7c)  may  be  rewritten  as 


|>,fl+lFc^=Ao=rfcr!S  (8) 
2  da,  £  2  c/a  ar 

Dividing  (8)  through  by  and  multiplying  by  r  to 
make  the  expression  nondimensional, 

1  P  „dG,,S,  ,1  ™dCmK_Jtt 

T2k?K~dZ;Sv'*+2tSE?~d^£ia-^dt‘  (9) 


Letting  ~jg  =  / 1  (reference  3) 


1  l2  dCZ{St 


1  Ic  dCm 


rfa,  6’  T,,?  +  2M  fcr* 


For  convenience  let 


J*  JT7/  1  /2  dCz.St 

™  (after  Glauert,  reference  3)  rj<7r  r^>K~, — 

Vr  df  2  A,>2  da,  6 

-  CLAy-2C„A  T'  — '^Aa=2rdU  «  “d  1  Jc_dC3_ 

2  ky1  da 

CcA7-2CtAF'-^A«=-2r^  b  1(7)  ,«»  ,,n 

rt0£  mffrfl-bM"l0Aa=T*  (11) 

lVc^i+\Vc<^Aa=Tkrin2  C  Since  9=(«+t)  and  (sec  fig.  1),  q  and  8  may 

Equations  (7a)  and  (7b)  are  nondimensional  and  tho  be  replaced  by  a  and  7  reducing  the  number  of  vari- 
onlv  variables  are  Ay,  AV and  Aa.  In  order  to  re-  ables  in  equations  (7a),  (7b),  and  (11)  to  three.  These 
duce  t7c)  to  a  comparable  form  it  is  necessary  to  .  equations  are  rewritten  as, 


rewrite  which  is  the  only  constant  in  the 

dq 

equations  not  in  readily  usable  form,  and  to  express 
6  in  terms  of  7  and  a. 

The  greater  portion  of  the  change  in  pitching  mo-  j2y\ 

ment  produced  by  rotation  occurs  because  the  rota-  dfi+di2)  c 

tion  changes  the  angle  of  attack  of  the  tail  surfaces.  ... 

A  positive  rotation  q  causes  the  tail  to  move  downward  The  foregoing  equations  of  equilibrium  must  each 
with  a  velocity,  relative  to  the  center  of  gravity,  equal  be  satisfied  at  each  instant  (neglecting  approximations 

to  ql  wffiere  l  is  the  distance  from  the  center  of  gravity  assumed)  of  flight.  It  is  obvious  that  each  of  the 

to  the  mean  quarter-chord  point  of  the  horizontal  variables  A7,  Aa,  and  AV  will  affect  each  of  the  others, 

surfaces.  The  tangent  of  the  change  in  angle  of  and  it  is  reasonable  to  assume  that,  if  any  one  of  them 

attack  at  the  tail  is  tan  A at=ql/V  which  may  be  re-  follows  a  regular  scheme  of  variation  with  time,  then 

placed  by  A a,=ql!V  for  the  small  angles  under  con-  the  others  will  vary  according  to  the  same  scheme. 


equations  arc  re  writ  ten  a: 


-CLAy-2CDAV'-^Aa=2i 


CDAy—2CLAV‘ 


-3*--* 


4 
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Assume  therefore  that  the  variables  change  accord¬ 
ing  to  the  exponential  expressions: 

A7  =  zlex< 

Aa=Bc*1 

±Y'=Ce*‘ 

(In  mathematical  terms  eu  is  an  integrating  factor.) 
where 

A,  B,  and  C  are  constants  depending  on  the  mag¬ 
nitude  of  the  initial  disturbance. 
e,  the  base  of  natural  logarithms  =  2.7 1828. 

X,  an  arbitrary  constant. 
t,  time  in  seconds. 

It  appears  that  if  X  is  positive  Ay,  Aa,  and  AU'  will 
increase  with  time  and  therefore  the  motion  will  he 
unstable.  If  X  is  negative  the  departures  from  the 
condition  of  steady  flight  will  decrease  and  the  motion 
will  be  stable. 

Now  if 

Ay=Aext 

etc. 

Substituting  in  (12) 

-CLAex'-2CDCeu-^Be’',=2T\Ceu  a 

Afi 

CDAeu-2CLCexl~~^B^tr=  -2  tA\pu  b 


m?r(BXex,4“  A\exc)-\- nmaBext  =  r(\2Bext-r  X2^4cx0r 
Each  of  the  equations  may  be  divided  by 
-  CLA-2C,>C-~^B=2,\C  a 

CDA-2CLC-dj£B=-lT\A  b 

mffrX(A+B)  (rX)2(^i  +  J5)  c 

From  (14b) 


(13) 


(14) 


2t\A+CdA-^-B 

da 

—'Tr, 


Also,  since  r  and  X  do  not  appear  except  in  the  product 
rX,  let  y  —  t\ 

Substituting  in  (14a) 


^B 


-cl(2  '■■a+C'A-Vqb) 

Multiplying  by  CL 

-  CcA- 2CDXA-  Cd-A + cJ^B  -  C'^B 

An 

=2(X')vi+  \'CdA-  y  ~[£B 


or,  since  Cl;-\-CT>2=CR1 

A[  -j-  CR2  4-  3  Cn  X'  4-  2  ( X' ) J] 

:+A£) 

Therefore, 


A 


rdCJ:_rdCp,ACL 


B  CJ+WoW +2{\'Y 
Also  from  (14c) 

A[m?X'  —  (X')2]  =  /?[(X')2— m?X'  —  nma] 
Therefore, 

A_tima-\-mg\'  —  (X7)2 
B  —  mtfX/  +  (X/)2 

It  therefore  appears  that  equilibrium  at  all  times  is 
possible  when, 

p  dCfj _ p  dCo  .  .  ,dCL 

D  da  yL  da  da  M^a-hn^X' —  (X')2  (15) 

ft2 +3  Co  X'  -j~  2  (X') 2  -m<X'-f(X')2 

or 

iV)r(c<g-ag>+x%) 

=  [C«2+3CroX'  +  2(X')2][Mma-hmtfX/  —  (X')2] 
Expansion  of  this  relationship  and  collection  of  terms 
give  the  biquadratic  expression 

(X' 


)4+(X')3[ 

[wB+l- 

xtl-™< 

i(s  cn+ 

ACC\_ 

da  ) 

(p  dCL 

A  Cd~^~ 

-C^+cAj\ 

r  i/ 

'  AC,, 

dC i) 

'  l"m^\ 

C°di~ 

C '-HI- 

—  tima~  =  0 


(16) 


For  simplicity  the  biquadratic  may  be  expressed  as: 
A(X,)4  +  5(X')3-f'C(X/)2H-/>X/-f£'=0  (17) 

where 

b-[-m*+K3Co+S)] 

DAprnyc^-cLd^+cA)-\cD^] 

p  2 

Unfortunately  there  is  no  simple,  direct  method  of 
solving  biquadratics.  It  is  possible,  however,  to  factor 
a  biquadratic  into  two  quadratics,-  each  of  which  is 
susceptible  to  direct  solution. 
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The  biquadratic 

(\T+£(x')J+a*')J+Dx'+/i'=<> 

may  be  replaced  by  the  expression 

[(X')s+a,X'  +  6,][(X')s+«jC  +  W  =  0 


from  which 


X'— 3±- 


It  appears  that 


/?=(ll  +  fl2 
<7=ai02  +  6l  +  62 
Z?=a!62+«2^i 
E =  6i62 


The  general  case  may  be  worked  out  as  follows: 
0.2 = B — 0 1 

C—  a,  (5— <ii)  +  61  +  62 

1^=0162+  — ®j)6i 

h-E 
62— r 


Therefore, 


C=ax{B— ad  +  6i  +  r- 


l?==Oi^-+ (5 — Oi)6i 
Dropping  the  subscripts: 


■W(D' 


VB-bD 


-C+6  +  - 


These  relations  can  be  solved  by  plotting  the  curves 
of  6  against  a.  There  are  two  intersections  of  these 
curves,  in  general,  corresponding  respectively  to  Oj, 
6t,  and  to  o3,  62. 

Also  from  (10) 

~=C-b-a(B-a )  a 

D—aC— ab—'a?{B— a)  +  6(fl— «)  b  (21) 

,  D-aC±a2B-(i' 
b~  B  —  'ln  ° 


By  use  of  the  foregoing  relationships  the  coefficients 
of  the  quadratics  may  be  determined  with  as  high  a 
degree  of  accuracy  as  desired  by  graphical  means  or 
by  trial  substitutions. 

At  first  glance  it  appears  that  a  (ith-power  equation 
such  as  (23)  would  be  harder  to  solve  than  a  4th-power 
equation  such  as  (17).  In  equation  (17),  however,  the 
complex  roots  must  be  obtained  ;  whereas  in  (23)  it  is 
necessary  to  solve  only  for  the  real  values.  Equation 
(23)  is  useful  from  a  practical  standpoint  chiefly  in 
obtaining  accurate  values  of  a-2  by  making  trial  sub¬ 
stitutions  fro  n  approximate  values  obtained  from  the 
expression  for  a,  given  on  page  6.  Because  of  the 
very  small  value  of  a2  it  is  generally  not  necessary  to 
include  the  terms  in  (23)  that  contain  powers  of  a2 
higher  than  the  third. 

Significance  of  X'. — As  appears  in  equations  (18), 
there  are  possible  either  4  real  values  of  X',  2  real 
and  1  pair  of  complex  values,  or  2  pairs  of  complex 
values.  The  values  of  B,  C,  and  D  in  the  normal 
flying  range  of  conventional  airplanes  are  always 
positive  because  of  the  signs  and  magnitudes  of  their 
constituent  factors.  It  is  obvious  that  no  positive  real 
value  of  X'  can  satisfy  the  biquadratic  unless  E  is 
negative  but  that  if  E  is  negative  there  is  such  a 
solution  for  X'.  A  positive  real  value  of  X'  signifies 
an  aperiodic  divergence.  If  —iima  is  positive  corre¬ 
sponding  to  static  stability,  then  E  is  positive  and  the 
biquadratic  expression  indicates  no  possibility  of  an 
aperiodic  divergence. 

The  values  of  /?,  C,  D,  and  E  are,  in  general,  such 
that  the  solution  for  X7  gives  two  pairs  of  complex 
values.  It  can  be  shown  by  mathematical  reasoning 
not  essential  to  this  treatment  that  an  expression  of 
Kex’r — where  X'  =  ±e'+  and  where  K,  j*' ,  and 

are  constants  and  i=y  —  1—  can  be  replaced  by  an 
equivalent  expression,  K' cos  (\f/'t'—6)  where  K' 
and  <5  are  new  constants.  Therefore  Ay  =  /1  eyr  may 
he  replaced  bv  Ay=A'et'r  cos  {ft’  —  5)  and  similarly 
for  the  other  variables.  It  appears  that  the  motion 
indicated  by  a  complex  root  is  therefore  made  up  of 
sinusoidal  variations  of  the  angle  of  attack,  angle  of 
altitude,  angle  of  the  flight  path,  and  velocity  along 
the  flight  path  and  that  the  amplitudes  of  the  oscilla¬ 
tions  increase  or  decrease  with  time  depending  upon 
whether  f'  is  positive  or  negative.  From  equations 
(18)  it  appears  that 


Substituting  (21)  in  (19a) 


H,  x  .  /Z)-aC+a23-rt3\ 

7=aCB_a)  - ) 

_l  jr( _ _ \ 


h'— — + 


h'  = 


o"— a5(3B)  +ft‘(3ffi+2C) 
—nHB,JriBO+(ii(2B!C+Bn+C--AE') 

-  a  (B'-D  4-  BV-  -  4  £B)  -f BCD  -D--B:E= 0  123) 
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where  f/,  jY,  and  \ p2'  express  the  stability  charac¬ 
teristics  of  the  airplane. 

From  the  expression  Ay  —  A'et'1'  cos  {ptf  —  b)  the 


period  of  the  oscillation  is  P' and  the  time  for 
the  oscillation  to  damp  to  one-half  amplitude  from  any 
instant  chosen  as  the  time  of  starting  is  T'  =  ~~  — 
Since  X'  was  set  equal  to  Xr  during  the  derivation,  Pf 
and  T'  are  in  units  of  time  which  are  equal  to 

seconds  or  the  time  in  seconds  necessary  for  the  air¬ 
plane  to  travel  the  length  l  at  the  velocity  V  multi¬ 
plied  by  the  relative-density  factor  y.  Therefore  the 
final  expression  for  the  time  in  seconds  for  the  oscilla¬ 
tion  to  decrease  to  one-half  amplitude  is 


-0.693  -0.693  m 

e  T~  v 


-0.313 

1  r' 


(24) 


under  standard  conditions;  and  the  period  in  seconds  is 


_2.83  \Wn 

-  r  \sCl 


(25) 


(BCD—D2—B2E)  is  negative.  This  fact  was  first 
pointed  out  by  Routh  (reference  4)  by  a  somewhat 
different  derivation  and  the  factor  is  known  as  “  Routh’s 
Discriminant.” 

There  has  been  developed  in  reference  L  an  approxi¬ 
mation  for  the  case  when  D  and  E  are  small  with 
respect  to  B  and  C  by  assuming  that 


ay—B 

bi~C 

Then 

D=Bb2+a.C 

E—Cb2 

or 

and 

D  BE 

a-i-Q  (jt 

It  follows  that, 


r.'=-§  (28) 


(29) 

(30) 


under  standard  conditions. 

Substituting  for  a  in  equation  (23): 


(31) 


(2nfl+ (2n*(3£) + (2n  w+2<7) 

-f-  (2p)'<JP+*BC)  4-  (2p)2(2B2CPBD+C2-4E) 

4-  (2p)  (B2D+BC2-4EB) +BCD- D*-B*E=Q 

(26) 

From  equation  (21) 

Reference  to  equation  (26)  reveals  that  from  the 
values  of  B,  C ,  D,  and  E  in  the  normal-flying  ranges  the 
coefficients  of  the  terms  in  2p  will  all  be  positive. 
Therefore  a  positive  real  value  of  p  can  exist  only  if 


These  approximate  expressions  are  generally  used 
and  arc  satisfactory  for  most  cases  encountered. 

Since  B  and  C  are  large  in  the  normal-flight  range, 
values  of  f,'  and  define  a  very  heavily  damped 
oscillation  which  is,  in  general,  of  low  period.  This 
oscillation  is  unimportant,  except  possibly  in  special 
cases  above  the  stall;  throughout  the  remainder  of  the 
report  p  and  will  be  used  without  the  subscripts 
to  refer  to  the  slightly  damped  phugoid  oscillation 
previously  defined  by  jV  and  fa'.. 

Derivation  of  expression  describing  the  sinusoidal 
motion. — The  expression  Ay^Ae^^Ae*'*’  can  be 
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replaced  by  the  equivalent  expression  A 7=^'"  where  A 
cos  (1 p't'—8)  where  A'  and  5  depend  upon  the  instant  variable 
from  whence  time  is  taken  as  zero.  If  time  is  zero  time  am 

when  Ay  is  at  a  point  of  maximum  amplitude  A'  =  Ay0  with  res 

and  5  =  0  so  that  These 

A7=Ato«'"' cos  ft’  (32)  tu‘ed  1. 

and  \p  . 

It  follows  that  f 

Aa=A<*0<d'''  cos  61)  (33) 

AVr/  =  AW<d'r  cos  ty't'—Bi)  (34)  checked 

Ad^Ad^'1'  cos  53)  (35)  and  tan 

/  ~  dCL  q  dCp  1  y.,dQ\  § 

Ayo„  da)Smd 

Aao  ^(3CD  +  4p) 

tan  6'=7ZT3K  „  dCD  ,  M,:\  &Ud±M 


where  Aa0,  AVV,  and  A 0O  are  the  magnitudes  that  the 
variables  would  have  were  they  at  a  maximum  at  zero 
time  and  5j,  82,  and  5^  are  phase  angles  of  the  variables 
with  respect  to  the  flight-path  angle. 

These  values  of  Ay,  Act ,  and  AVf  may  be  substi¬ 
tuted  in  equations  (12)  and  solutions  made  for 
and  \p'.  It  will  be  found  that  the  complete  solution 
for  £'  and  given  in  the  preceding  section  can  be 

checked  in  this  manner  and  the  values  of  tan5,» 

and  tan  6,  are  found  at  intermediate  steps  to  be: 


d(\ 

+(3c^4ncosa' 


'(*#- 


cze+rB)-*?, 


{mA'+m*- Kf')2— fiA')2]} 


n  I  y/dCp\  dCi , 

1E7  -  a  da 

—m, - sm  cos  5., 


2fCLd-z£{2!’-m,)- 


_ HU  _  _ 

-4r'm,-2(^'),t-f)(n!!  +  (2f'-/nt)(cu-|-^)] 


The  detailed  derivation  is  not  repeated  because  it 
is  quite  long  and  tedious  and  introduces  no  new  con¬ 
cepts. 

Application  of  mathematical  formulas— The  rela¬ 
tionships  derived  in  the  preceding  section  of  this 
report  make  possible  various  treatments  of  the  prob¬ 
lem  of  determining  stability.  The  period  and  time 
to  damp  to  one-half  amplitude  may  be  determined 

■  /t  dCfl  dCp  i 

directly  having  given  CLi  G0>  ana 

—  Mma  from  the  equations  (17),  (24),  (25),  and  either 

(20)  and  (27)  or  (30)  and  (31). 

An  understanding  of  the  underlying  principles  gov¬ 
erning  stability  can  be  had  only  from  a  consideration 
of  the  variations  of  the  angle  of  attack,  angle  of  atti¬ 
tude,  angle  of  flight  path,  and  velocity  along  the  flight 


path  relative  to  each  other.  Figures  2  to  14  are  in¬ 
cluded  to  show  the  nature  of  these  variations.  Figure.; 
2  to  9  present  the  effects  of  changing  —  mq  and  —  fxma 
independently  for  two  sets  of  typical  values  of 

CL}  CDl  and  Figures  10  to  14  illustrate 

qualitatively  the  variations  in  phase  relationships  re¬ 
sulting  from  changes  in  — ixtna  and  serve  as  a  basis 
for  the  discussion  of  their  effects  upon  the  stability 
of  the  airplane.  The  relationships  (36)  to  (39)  can 
be  used  to  determine  the  characteristics  of  the  motion 
for  any  particular  design. 

The  solution  lor  the  stability  characteristics  gives 
no  direct  indication  of  the  effect  of  variation  of  the 
individual  parameters  upon  the  stability  characteristics. 
By  means  of  the  relationship  that  follows  directly 
from  equation  (22) 


Tnq{cJr2^)Jrnma-~[d+2^cJt^(^)2]  k 

,  (-mg[^+2r/g+4(r)2l-iumg(e+2r/)+2r/^+2r/c+4(r/)2]l 
+  i  —  wa-b(c+ 4f')  ) 

{  J _ _ — ma+(c+ 4f') _ 1 

>m«|^m?[(/+2^c+4(n2l-^(6+2r,)+2r,y+2r'c-h4(r')2ir 


121002-35 - 2 


8 

where 
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charts  may  be  plotted  showing  the  variations  of  —  mq 
with  — uma.  necessary  to  secure  given  values  of 
and  \p' .  A  number  of  these  charts  have  been  prepared 
covering  the  range  of  conditions  likely  to  be  encoun¬ 
tered  in  normal  flight  and  are  included  in  this  report. 
(See  figs.  15  to  54.) 

DISCUSSION 

The  mathematical  relationships  evolved  in  the  pre¬ 
ceding  paragraphs  permit  calculation  of  the  probable 
stability  characteristics  of  a  proposed  design,  but  they 
offer  little  information  as  to  the  relative  importance 
of  various  factors  or  as  to  the  reasons  for  the  effects 
produced  by  changes  in  those  factors.  In  the  following 
paragraphs  the  oscillatory  motion  is  first  considered 
in  detail  so  that  the  sequence  of  flight  conditions  that 
must  exist  if  instability  is  to  arise  may  be  pointed  out. 
Next  is  given  a  general  discussion  of  charts  (figs.  15 
to  54)  that  show  the  effects  on  the  stability  character- 

dC 

istics  of  the  six  fundamental  parameters  CL,  Co, 

—uma,  and  —to..  Finally,  the  effects  of  various 
da 

physical  characteristics  of  the  airplane  on  its  stability 
are  considered  in  the  light  of  the  earlier  discussion. 

The  oscillatory  motion.-— When  the  terms  in  Al" 
were  eliminated  by  simultaneous  solution  of  the  equa¬ 
tions  of  equilibrium  tangent  and  normal  to  the  flight 
path,  respectively,  the  following  equality  was  found  to 
exist: 


p  dCL _ p  dCp  .  ydCr 

D  da  L  da  *  da  groa+m^X'  — -  (X')2 

'  C r  -r  3 CD\f  -f  2  (X ' ) ■  " [\'y 

The  terms  to  the  left  of  the  equality  sign  arise  out 
the  necessity  for  equilibrium  of  forces;  the  terms  to  t 
right  arise  out  of  the  necessity  for  equilibrium 
moments.  If  a  case  be  investigated  with  —  y.ma  —  ( 

C^-CL^+yd-^=-CR*-ZCDW-2(X)-‘ 


or 


Since  both  CD  and  are  positive  below  the  stall. 

follows  that  if  — is  zero  the  linear  motion  of  t 
airplane  is  a  heavily  damped  oscillation  with  the  peri< 


P'  =  - 


2t 


v,I(r-S-'!‘1&+';-,)-ro(3a+©' 

and  the  time  to  damp  to  one-half  amplitude 


T 


0.693 

K3C"+© 


all  in  nondiinensional  units,  or  the  motion  is  an  apei 
odic  convergence  with  the  time  to  converge  to  one-lm 
amplitude 


0.693 


-o^.+cA 


depending  on  the  relative  magnitudes  of  the  quantities. 

Since  instability  of  the  linear  motion  cannot  arise 
when  —pm*  is  zero,  it  seems  probable  that  under¬ 
standing  as  to  the  underlying  causes  of  instability  may 
be  gained  by  considering  the  effect  of  — /rm«  upon 
the  sequence  of  flight  conditions  during  oscillatory 
motion.  In  the  following  paragraphs  the  part  played 
by  —ixTYia  is  analyzed  by  physical  reasoning.  Only 
the  case  when  —  ixma  is  positive  is  considered  because 
when  —  fxma  is  negative  an  aperiodic  divergence  from 
the  equilibrium  condition  occurs  and  the  question  of 


dynamic  stability  does  not  arise,  it  is  iirst  pointi 
out  that  —  fima  plays  a  primary  part  in  detenninii; 
the  phase  angles  between  the  change  of  angle  of  attac 
and  the  changes  in  attitude  and  flight-path  angle 
It  is  next  shown  that  stability  of  the  angular  inotio 
depends  upon  the  phase  angle  between  the  angle-o 
attack  change  and  the  attitude  change.  It  is  finall 
shown  that  instability  of  the  linear  motion  can  arb 
only  when  the  phase  angle  between  the  angle-o f-attac 
change  and  the  flight-path-angle  chauge  falls  within 
certain  range. 
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Figures  10  to  12  illustrate  qualitatively  the  effect 
of  —  ixma  upon  5,  and  8if  the  phase  angles  between  the 
angle  of  attack  and  the  flight  path  and  attitude  angles, 
respectively.  When  —p ma  is  zero  the  attitude  does 
not  change  with  time  and  hence,  if  an  oscillation  is  set 
up  in  the  flight-path  angle  (7),  it  is  obvious  that  the 
angle  of  attack  must  change  as  shown  in  figure  10  where 
5,  is  180°.  Although  not  shown  in  figure  10,  physical 
reasoning  leads  to  the  conclusion  that  64  is  90°.  When 

—  pma  is  very  large  the  moment  tending  to  cause  the 
attitude  to  change  with  the  flight  path  is  also  large 
and  the  angle-of-attack  change  is  small.  The  condi¬ 
tions  are  approximated  by  figure  12.  The  phase 
angles  Si  and  54  are  also  small  and  approach  zero  as 

—  fjuna  is  increased.  Figure  11  shows  an  intermediate 
case  which  more  truly  represents  the  usual  condition 
than  either  of  the  others.  Here  5|<^180°  and  54<90°. 
The  curves  of  figures  10  to  12  were  drawn  from  physical 
reasoning  but  are  fully  confirmed  by  the  computed 
curves  of  figures  2  to  4  and  7  and  8. 

Referring  to  the  curves  of  figures  2  to  9  one  sees 
that,  if  2  points  are  chosen  exactly  1  period  apart  when 
A7  is  zero,  the  value  of  A  V'  will  be  found  to  be  greater 
at  the  second  point  than  at  the  first  if  the  motion 
is  unstable,  to  be  unchanged  if  the  motion  is  neutrally 
stable,  and  to  be  less  if  the  motion  is  stable.  Similar 
observations  can  be  made  with  respect  to  the  val¬ 
ues  of  Ay  when  AV'  is  zero  and  to  the  values  of 


^=-when  Ad  is  zero.  In  other  words,  when  neutral 
at 

stability  exists,  the  momentum  along  the  average 
flight  path  m(V+  A 10,  the  momentum  perpendicular 
to  the  average  flight  path  {mV sin  A7),  and  the  angular 
momentum  are  each  the  same  at  the  end  of  a  period 
as  at  the  beginning.  If  instability  exists,  there  is 
more  momentum  at  the  end  of  the  period  for  each  of 
the  types  of  motion  than  at  the  beginning.  It  is 
therefore  desirable  to  consider  the  sequence  of  flight 
conditions  that  results  in  the  increase  or  decrease  of 
momentum  in  the  system. 

Consider  first  the  angular  motion.  Figure  13  is  an 
illustration  of  a  typical  set  of  variations  of  Aa,  A7, 
and  Ad  when  f'  is  zero.  At  the  beginning  of  the  half 


period  when  A0  is  a  maximum,  A q  (^or  A  is  zero  and 


the  angular  momentum  BAq  is  therefore  zero.  At 
the  end  of  the  half  period  Aq  is  again  zero.  It  there¬ 
fore  follows  that  the  angular  momentum  introduced 
into  the  system  during  the  half  cycle  must  be  balanced 
by  the  momentum  removed  during  the  same  interval. 

At  each  instant  when  there  is  a  finite  value  of  Aa, 
there  is  a  pitching  moment  introduced  into  the  system 


equal  to  Aa\pV2Sc^^ 


Since  a  moment  M  applied 


during  a  time  t  introduces  angular  momentum  into  the 
system  equal  to  Mt 


mk 


<sw 


1W-0 

fit=w  1  d( 7 

=  (  —  pV2Sc—r-^ Aa0ef*  cos  (1 pt- 

J<H= 0  -  aa 


a ApV’Sc^dt 


da 


■St)dt 


represents  the  angular  momentum  introduced  into  the 
system  during  one-half  cycle  (neglecting  the  effects  of 
terms  involving  the  products  of  two  small  quantities, 
as  has  been  done  throughout  the  discussion). 

Considering,  for  simplicity,  the  case  of  neutral 
stabilitv 


mkKft)rSl  lhpv'8eife*a' cos 

T" 


1  T/2  O  A  dC„ 

2fiVSc 


sin  {\f/t— $4) 


pVs&Aa„^fsm  a. 


is  the  angular  momentum  introduced  into  the  system. 

The  rotational  velocity  causes  a  change  of  angle  of 
attack  at  the  tail  that  introduces  a  moment  tending  to 
oppose  the  rotation.  The  magnitude  of  this  change  is 
,  d&  -,dd 
=l**=  di=  di 

since 

.dd_dd 


The  moment  is  —  Kr)t^pV2S  r^- 


1 


..de 
S,\ldt 
sJ  V  • 


(See  p.  3.) 

The  angular  momentum  removed  from  the  system 
(f=0)  during  the  half  cycle  is 


mk} 


Tn  the  case  of  neutral  stability  the  momentum  intro¬ 
duced  must  equal  that  removed  and  therefore 
dCm 

sin  «,  /<!(’?, S,  s. 

i  -’MVpl'.S :l\daix)^u 

or 

pma  Aa0 

V  A0oS 

From  figures  2,  3,  4,  7,  and  8  it  is  evident  that  as 

—p.ma  increases  i-j  and  64  all  decrease.  The  first 

Y  A70 

term  in  equation  (41)  may  therefore  either  increase 
or  decrease  with  increase  in  —pina‘,  hence  increasing 
—  ju ma  may  decrease  the  angular  momentum  added 
during  a  cycle.  What  happens  in  an  actual  case 
depends  on  the  magnitudes  of  the  quantities,  as  wdll  be 
brought  out  later  in  the  discussion. 


T7  -r^-sin  m,=0 


(41) 
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Figure  14. — Variations  of  angle  of  attack,  rate 
angle  of  flight  path  for  an  illustrative  case. 


of  change  of  angle  of  flight  path,  ami 
Damping  coefficient  assumed  zero. 
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The  introduction  of  linear  momentum  along  and 
perpendicular  to  the  flight  path  might  he  discussed  by 
considering  equations  (l 2a)  and  (12b)  separately  but 
the  analysis  may  be  simplified  by  combining  the  two 

through  elimination  of  terms  involving  AV'  and  — r~- 

at 

The  relationship  established  is: 


thought  that  the  effect  of 


~r~~~  as  an  unstabilizin 
da 


-CJly- 3C„r'£— 2  AX- 


fy  dCj,  S~1  dC O 
C°da  tLda 


dCL  da 
7  da  d  t 


Considering  the  case  of  zero  damping  during  the  in¬ 
terval  from  ^f  =  0  to  \pt~T  it  is  obvious  from  figure  14 
that  the  integral  sum  of  terms  multiplied  by  A7  or  by 

— /  is  zero  because  these  quantities  are  symmetrical 

with  respect  to  the  zero  line.  The  equation  of  equality 
of  momentum  introduced  to  that  removed  becomes: 


p--»&  ru. 

J  *t= 0  dt  j  0  V 

rm\M-  d«dt 

J  a- 0  da  dt 


c^-°^)Aadt 


which  reduces  to 


-3  Co- 


Aao  •  .  .  dCL  Aa0  *  n 
■7~~  sin  61  ^  cos  5l=0 

A7q  da  A70 


The  growth  of  linear  momentum  therefore  depends 
upon  <5p  and  \p'  in  addition  to  the  parameters 

CDj  CL>  and  Since  is  normally  positive 

the  following  facts  appear:  CD  is  a  factor  serving 
always  to  remove  momentum  from  the  system; 

-  is  a  factor  adding  energy  to  the 
system  when  as  is  normally  the  case 

at  high  angles  of  attack;  and  adds  momentum  to 

the  system  when  5,  <90°  and  removes  momentum 

when  5('  0(1°.  When  ii/na  is  zero,  or-  IN0°  and  tin*  j 

factor  removing  energy  from  the  system  is  ^36’/,+^^  ^ 
as  was  shown  earlier  in  the  discussion.  As  —y.ma  is 

increased,  and  ~  all  decrease  (figs.  2.  3,  4.  7. 

A  7,]  y 

and  8).  The  factor  sin  ,, 

reaches  its  maximum  at  some  point  between  5!  =  1 80° 

and  <5t  =  90°.  The  factor  cos  <5t  changes  from 

da  A7o 

a  stabilizing  to  an  unstabilizing  influence  at  5,  =  90° 
and  reaches  a  maximum  as  an  unstabilizing  element 
between  S^OO0  and  5i  =  0.  It  would  appear  at  first 


factor  would  always  be  small  because  of  the  decrease  <> 

-r^with  increase  of  cos  <$!  but  such  is  not  necessarib 
A70 

the  case  because  of  the  large  ratio  of  to  CDf  particu 

larlv  at  low  angles  of  attack. 

There  is  an  interesting  point  in  connection  with  th 

factor  —  that  appears  as  the  ohie 

unstabilizing  element  at  small  values  of  If  th 

law  of  induced  drag  he  applied  to  determine  Cn  am 

dCD 

it  appears  that 

(c  dCL  pdCn\_r  dCl.p  dCL  p  dCa, 

V  “  ’■  Ta  )-' (  °’d^ +Cn'da  Cl1U 


_p  dCL,K 
~Cl>Pd^+- 


_dCL(r  CL\ 

t ) 

Therefore  equation  (42)  may  be  rewritten  as 

~3°n+^  ^^Cn-Gnp)j,  sin«,+cos5,J=0  (4:t 

It  appears  that  instability  of  the  linear  motion  foi 
small  values  of  —iitna  can  occur  only  at  angles  ol 
attack  where  CDi^>CDp.  The  detrimental  effect  of 
Cn{  is  not  as  great  as  at  first  appears,  however,  be¬ 
cause  the  damping  factor  3  CD  increases  with  Cn{  also. 

In  the  foregoing  analysis,  the  effect  of  f'  has  been 
neglected  and  therefore  the  statements  made  camim 
be  considered  rigorously  true  when  the  stability  or  in¬ 
stability  is  of  appreciable  magnitude.  The  analysis 
does,  however,  point  out  the  more  important  influences 
and  the  nature  of  the  interaction  of  events  that  brings 
them  into  play.  Equation  (43)  cannot  he  applied  near 

26’/— 

the  stall  where  ceases  to  equal - g^-- 

In  figures  2  to  14,  inclusive,  and  in  the  discussion  of 
the  oscillatory  motion  an  arbitrary  value  of  A70=  i  has 
been  assumed.  No  attempt  will  be  made  to  evaluate 
Ay0  in  terms  of  control  movements  or  gust  velocities 
as  the  actual  magnitude  is  unimportant  with  respect 
to  stability  (assuming  that  the  deviations  are  not  so 
great  as  to  destroy  the  validity  of  the  basic  assump¬ 
tions).  It  seems  not  out  of  place,  however,  to  suggest 

that  relative  magnitudes  of  5l}  62,  etc., 

^7o  A7o  A70 
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may  have  important  bearing  upon  the  comfort,  ease  of 
handling,  and  design  load  factor  of  aircraft  and  that 
further  work  taking  these  factors  into  consideration 
might  lead  to  valuable  information. 

Charts  of  rotational  damping  factor  against  static 
stability  factor.— In  the  section  dealing  with  applica¬ 
tion  of  mathematical  formulas  a  convenient  graphical 
means  of  showing  the  variation  of  and  \J/'  with  —mq 
and  —pma  lias  been  described.  This  method  of  pres¬ 
entation  was  given  by  Gates  (reference  5)  in  essen¬ 
tially  the  same  fashion  although  differing  in  detail. 
In  reference  5  preference  is  given  to  charts  with 
coordinates  of  tail  volume  and  fore-and-aft  location 
of  the  center  of  gravity.  Preference  has  been  given 
in  this  report  to  charts  with  coordinates  of  rotational 
damping  factor  and  static-stability  factor  because  they 
are  more  convenient  for  use  with  data  from  wind- 
tunnel  tests,  they  do  not  require  assumption  of  arbi¬ 
trary  fixed  values  for  several  important  factors,  and 
they  cover  a  much  wider  range  for  a  given  number  of 
charts.  The  charts  are  entirely  nondimensional  and 

vary  only  with  CL ,  CDi  ^ j and  \ariation  of 


,  dot  i  , 

empirical  factors  such  as  the  K  m  vu  ^  \see 

p.  17),  etc.,  affect  only  — mtf  or  —  as  the  case  may 
be,  and  their  effects  upon  the  stability  are  readily 
apparent. 

The  series  of  charts  presented  (figs.  15  to  54)  are 
intended  to  show  somewhat  more  precisely  and  com¬ 
pletely  than  has  been  done  in  the  preceding  discussion 

„  „  dCt  ■  A  dCn 

the  effects  of  the  parameters  CL}  CD,  and 

upon  dynamic  stability  and  to  provide  a  convenient 
graphical  means  by  which  the  designer  may  estimate 
the  probable  stability  characteristics  of  a  proposed 
airplane  and  the  effects  of  various  changes  without 
recourse  to  extensive  calculations.  The  charts  cover 
the  range  of  values  of  the  parameters  that  appear 
likely  to  be  attained  in  the  near  future.  The  values 
represented  are  summarized  in  table  I. 

Certain  general  characteristics  of  the  charts  are  im¬ 
mediately  apparentupon  inspection.  As  —  n>na  increases 
from  zero,  f'  at  first  becomes  more  positive  correspond¬ 
ing  to  a  decrease  in  dynamic  stability  but,  at  a  fairly 
small  value  of  -  changes  its  trend  and  becomes  more 
negative.  This  tendency  is  general  throughout  all  the 
charts  and  is  apparent  whether  —mQ  is  large  or  small. 
It  will  be  remembered  that  such  an  effect  appeared 
probable  from  the  discussion  of  the  introduction  of 
angular  momentum  into  the  system  during  a  cycle 
K v  —m.  As  nointed  out  at  that  stage  of  the  report^ 


— ,  sin  5<t  and  all  decrease  with  increase  of  —^ma; 
A7o  v  .  . 

the  charts  of  variation  of  the  stability  characteristics 
with  —rnq  and  —pm«  show  definitely  that  the  decrease 
of  the  product  of  these  factors  is  much  more  than 
sufficient  to  nullify  the  increases  in  —  after  a 
certain  value  of  —  j ama  has  been  exceeded.  Increase 

121002—  ;jr> - ;# 


of  xf/',  hence  decrease  in  period  of  the  oscillations,  with 
increase  of  — /uma  occurs  at  practically  all  values  of 
—mQ  and  —  /im8. 

Increases  of  — mq  give,  in  general,  more  negative 
values  of  f '  and  consequently  more  rapid  dying  out  of 
the  oscillations.  At  large  values  of  —  iima  increases 
in  —mQ  give  more  positive  values  of  but  this  effect 
is  not  of  practical  importance.  It  will  be  noticed  that 
at  fairly  large  values  of  —  increasing  —  mq  has  hut 
slight  effect  but  that  — becomes  of  increasing  im¬ 
portance  as  -/iffla  is  made  smaller.  Increasing  —  mq 
decreases  V  rather  gradually,  giving  oscillations  of 
longer  period.  It  appears  that  if  the  criterion  of  sta¬ 
bility  be  taken  as  the  number  of  oscillations  necessary 
for  the  amplitude  to  decrease  to  one-half  its  original 
value  then  the  value  of  —  mq  is  of  particular  importance. 

The  charts  of  figures  23,  27,  28,  and  29  show  that 
increasing  Ch  without  changing  other  factors  increases 
quite  markedly  the  tendency  to  instability.  This  ob¬ 
servation  agrees  well  wdth  the  effect  to  be  expected 
from  increasing  CL  that  might  have  been  predicted  from 

the  part  played  by  CL  in  the  factor  ) 


in  equation  (42). 

The  effect  of  increasing  CD  without  changing  other 
factors  appears  in  figures  29,  30,  and  31.  It  will  be 
seen  that  increasing  CD  reduces  the  range  of  instability 
and  brings  the  curves  of  nearer  together  with  a 
very  large  net  increase  in  damping.  The  effect  upon 
\ff'  is  negligible. 

The  slope  of  the  lift  curve  is  of  importance  although 
not  giving  such  extreme  effects  as  changes  in  CL  and  CD 
(figs.  34,  35,  and  36).  Here  again  the  effect  is  in  good 

agreement  with  equation  (42).  Increasing  ~~  tends 


to  extend  the  instability  region  to  greater  values  of 
or,  in  other  words,  if  —  pm*  is  large,  is  less  than 

90°  an(i  becomes  a  definite  unstabilizing  factor. 
da 

At  small  values  of  -m„,  ^  is  a  stabilizing  factor, 
as  is  apparent  from  the  rurves.  Increasing  ^  also 


tends  to  decrease 

Increasing  the  slope  of  the  drag  curve  is  distinctly 
unfavorable  to  stability  as  shown  in  figures  42,  44,  and 
45.  The  period  of  the  oscillation  is  but  slightly 
affected. 

Reference  to  charts  15,  21,  and  33  shows  that  the 
general  effect  of  increasing  angle  of  attack  is  to  in¬ 
crease  the  region  of  instability  but  to  decrease  tin- 
spread  between  curves  of  £'=k.  It  appears  that 
instability  of  the  oscillatory  motion  is  very  unlikely 
for  smalf  angles  of  attack  and  becomes  increasingly 
likely  with  increase  of  that  angle.  On  the  other  hand, 
if  —^ma  and  —  mq  remain  unchanged  throughout  the 
change  of  angle  of  attack,  the  value  of  £'  may  become 
more  negative,  which  seems  to  be  directly  contra¬ 
dictory  to  the  preceding  statement.  The  explanation 
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lies  in  the  fact  that  at  small  angles  of  attack  changes 
of  —pma  and  —mq  have  little  effect  upon  the  damping 
coefficient;  whereas  at  large  angles  of  attack  a  small 
change  in  —  pma  or  —niq  may  change  the  damping 
coefficient  from  a  definite  negative  value  to  a  definite 
positive  value.  The  period  of  the  oscillation  decreases 
markedly  with  the  increase  of  angle  of  attack  but  the 
decrease  is  not  so  great  as  would  appear  to  be  the  case 
from  consideration  of  the  increase  in  i>' .  The  period 


■Jc 

in  seconds  is  proportional  to  —7 -  and  the  increase  in 


Ch  is  sufficiently  great  partly  to  counteract  the  increase 
in  yf/f . 

The  effects  of  various  physical  characteristics  of 
the  airplane  on  the  stability. — In  the  preceding  dis¬ 
cussion  no  consideration  has  been  given  to  the  various 
dimensional  characteristics,  aerodynamic  interferences, 
etc.,  that  determine  the  values  of  the  fundamental 
parameters.  A  large  number  of  factors  affect  the 
stability  but  in  many  cases  the  effects  are  of  a  minor 
nature.  Only  the  more  important  ones  will  be  dis¬ 
cussed  in  the  following  paragraphs. 

The  wing  loading  appears  in  but  two  places  in  the 
analysis  of  stability  in  this  report.  The  nondimen¬ 
sion  al  relative-density  factor 


m 

w 

s 


32.2  pi 

■'CD 
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(standard  conditions  at  sea  level) 


is  directly  proportional  to  wing  loading  and  the  factor 


m 

r~PSV~ 


2.22 


is  proportional  to  the  square  root  of 


wing  loading.  Since  m  appears  only  as  the  coefficient 

of  —ma  it  is  apparent  that  the  variation  of  the  non- 
dimensional  values  of  f'  and  i>'  with  g  are  the  same  as 
for  variations  of  —pma.  It  therefore  appears  that, 
if  the  longitudinal  motion  of  an  airplane  with  respect 
to  the  air  is  to  be  unaffected,  an  increase  in  the  wing 
loading  of  a  given  design  must  be  accompanied  bv  a 
dC 

proportional  decrease  in  provided  that  f,  c,  and 

kY  are  unchanged.  Aside  from  the  effect  upon  —pma 
the  wing  loading  affects  the  stability  characteristics 
through  the  time  factor  r,  which  appears  in  the  reduc¬ 
tion  of  the  nondimensional  units  $*'  and  to  the  time 
to  damp  to  one-half  amplitude  and  to  the  period  of 
the  oscillation  in  seconds 


and 


„  -0.093 


w 


If  —pma  is  preserved  unchanged  as  -g-  is  increased 


both  the  period  and  the  time  to  damp  are  increased 
but  the  number  of  oscillations  for  a  given  degree  of 
damping  is  unchanged.  The  effect  of  an  increase  of 
wing  loading  without  change  of  other  factors  is  there¬ 
fore  to  increase  both  ^  and  r  with  the  net  effect,  in 
general,  of  increasing  the  period;  to  make  f'  more 
negative  if  —pma  is  large,  with  a  net  effect  upon  the 
time  to  damp  either  of  an  increase  or  of  a  decrease; 
and  to  make  £'  less  negative  if  —pma  is  small,  resulting 
in  a  comparatively  large  net  increase  in  the  time  to 
damp  to  one-half  amplitude. 

The  aspect  ratio  of  the  wing  combination  affects  the 
stability  indirectly  through  its  effect  upon  the  para- 
rfft  dCD 

da  da 


meters 


and  CD  at  a  given  value  of  CL. 


Increasing  the  aspect  ratio  increases  and  decreases 


both  and  CD . 


Comparison  of  the  curves  of  figures 


36  and  39  corresponding  to  aspect  ratios  of  5  and  8, 
respectively,  shows  that  (neglecting  the  effect  due  to 


dC 

change  of  which  is  comparatively  small  and  will 

not  materially  affect  the  conclusions  to  be  drawn)  the 
net  effect  of  increasing  the  aspect  ratio  is  to  increase 


the  range  of  instability.  The  decrease  in  an 

effect  favorable  to  stability,  is  less  important  than  the 
decrease  in  CD)  which  quite  markedly  tends  to  decrease 
the  stability. 

The  effect  of  parasite  drag  is  brought  out  quite 
clearly  in  figures  33,  36,  and  40.  (In  general,  increas¬ 
ing  the  parasite  drag  also  increases  the  slope  of  the 
drag  curve  and  this  fact  has  been  taken  into  account 
in  the  figures  by  using  the  method  of  reference  6.) 
It  is  apparent  that  the  drag  is  an  important  item  in 
determining  the  stability  characteristics  of  a  design 
and  that  the  “cleaner'’  the  airplane  the  greater  the 
tendency  to  unstable  oscillations.  Increasing  the 
parasite  drag  results  in  a  decrease  in  the  unstabilizing 
factor  {CD{— CDp)  and  an  increase  in  the  stabilizing 
factor  3 CD,  which  appears  in  equation  (42)  in  the 
discussion  of  the  growth  of  linear  momentum  during  a 
cycle.  It  appears  therefore  that  air  brakes  of  various 
kinds  may  be  expected  to  have  beneficial  effects  upon 
the  longitudinal  stability  provided  that  they  do  not 
introduce  undesirable  moments  or  interferences. 


Highly  efficient  designs  that  achieve  high  values  of 
CLmax  without  the  use  of  devices  introducing  parasite 
drag  may  be  expected  to  be  deficient  in  longitudinal 
damping. 

The  moment  of  inertia  of  the  airplane  about  the 
lateral  axis  appears  as  the  term  kY 2  in  the  denominators 
of  the  expressions  for  ~mq  and  -ma.  Increasing  kY 2 
decreases  both  and  -pma  proportionally  with  a 

corresponding  increase  in  the  time  to  damp.  The 
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effect  may  be  readily  visualized  by  imagining  any 
point  —jj.}naf  —wiq  upon  the  chart  moved  from  or 
toward  the  origin  along  a  radial  line  an  amount  such 
that  the  ratio  of  the  new  distance  to  the  origin  to  the 
original  distance  is  the  ratio  of  the  original  value  of 
ky 2  to  the  new  value.  It  is  therefore  to  be  expected 
that  distribution  of  the  mass  of  an  airplane  along  its 
longitudinal  axis  may  lead  to  unstable  oscillations 
although  other  factors  such  as  length  and  size  of  tail, 
location  of  the  center  of  gravity,  etc.,  are  such  that 
ordinarily  the  damping  would  be  satisfactory.  The 
effect  upon  the  period  of  changing  JcY2  is  not  great 
although  increasing  kY 2  will,  in  general,  lead  to  oscil¬ 
lations  of  shorter  period. 

The  location  of  the  center  of  gravity  affects  the 
longitudinal  stability  through  its  effect  upon  the 

term ^2  appearing  in  the  expression  for  —  ma.  The 

pitching  moment  of  an  airplane  about  the  center  of 
gravity  may  be  expressed  as 


longitudinal  center-of-gravitv  location  upon  the  sta¬ 
bility  characteristics  thus  follows  from  the  preceding 
discussion  of  the  effect  of  changes  in  —nma.  The 
location  of  the  aerodynamic  center  of  the  wing,  de¬ 
fined  by  Cat  affects  the  stability  in  the  same  manner 
as  the  location  of  the  center  of  gravity  but  with  the 
opposite  sign.  The  value  of  Ca  given  by  wing  theory 
is  0.25  (reference  8).  Wind-tunnel  data  indicate  the 
actual  values  to  range  from  0.23  to  0.25  for  conven¬ 
tional  airfoils  (reference  9). 

The  vertical  location  of  the  center  of  gravity  with 
respect  to  the  mean  aerodynamic  wing  chord  is  not 
unimportant  in  many  cases  and  should  be  considered 
particularly  in  high- wing  or  low- wing  monoplanes. 

Since  -  is  positive,  positive  values  of  —  (as  in  the 

case  of  a  low-wing  monoplane)  make  less  negative 

and  may  cause  an  unstable  divergence  (static  insta¬ 
bility)  in  cases  that  would  be  considered  stable  on  the 
basis  of  calculations  involving  only  Cg,  The  value  of 


-CzJP'-Cd  +  Cm^+Cx^  +Cnp  ~  is  derived  as  follows: 


dCZt  l  St 

+“'■3^7  e  S”‘ 

Differentiating  with  respect  to  angle  of  attack 

<n  -r  )+d^?-+d^ 

da  da  0  a  da  C  da 

,  datdCztlSt 

'  da  dat  CS" 


,  cos  a-\-CL  sin  a 


cos  “+(°D»+©sin  “ 

where  Cdw  and  are  for  the  wing  alone.  For  a  case 


with  ct=1.4,^ 


=4,  CDo  =0.01,4=6,  a=15° 


or,  for  practical  purposes, 


*j»-gcc,-a>+ 


-a>+% 

datdCZt  l  St 


Us  IX  f l 

'  ^ 1  da  dat  c  S 

where 

Cg,  ratio  of  distance  of  c.  g.  from  the  leading  edge 
of  the  wing  to  the  wing  chord. 

Cai  ratio  of  the  distance  of  the  aerodynamic  center  of 
the  wing  from  the  leading  edge  of  the  wing  to 
the  wing  chord. 

CZirt  normal-force  coefficient  of  wing. 

CXtt ,  longitudinal-force  coefficient  of  wing. 

2,  distance  of  mean  wing  chord  below  c.g. 

C„  parasite  pitching-moment  coefficient  of  fuselage 

and  binding  gear. 

— ,  ratio  of  change  of  angle  of  attack  at  the  tail  to 
^ a  change  of  angle  of  attack  at  the  wdng,  an  em¬ 

pirical  factor  depending  upon  downwash  (refer¬ 
ence  7). 

since  is  positive,  increasing  C0  by  moving  the 

d<X  dCm 

center  of  gravity  toward  the  rear  tends  to  make 
W  npffativp  nnH  to  rleerense  —  u  .  The  effect  of 


Co„=0.01  +  gj-0.11 

dCD„  2X1.4X4_„,„ 


-4^=  (0.8X0.960)  +  (4.11  X0.259) 

da 

0.77-}-  1.06  =  1.83 

Assuming  the  case  of  a  low- wing  monoplane  with  -  — 
0.25  and  C0  =  0.30,  it  will  be  found  that  the  stability 
will  be  the  same  as  though 

?=0  and  C,  =  0.30+(^X0.25)  =0.41. 

The  example  represents  perhaps  an  extreme  case  but 
illustrates  the  desirability  of  investigating  the  effect  of 
vertical  center-of-gravity  location. 

The  value  of  the  slope  of  the  curve  of  parasite 
pitching-moment  coefficient  against  angle  of  attack 
may  bo  determined  with  reasonable  accuracy  only  by 
careful  testing  of  a  scale  model  or  the  complete  air¬ 
plane.  For  a  number  of  designs  of  military  airplanes 
for  which  sample  computations  were  made  and  com¬ 
pared  with  measured  values  of  values  of 

were  found  to  range  from  0  to  0.4,  with  an  average 
value  of  0.2,  at  angles  of  attack  corresponding  to 

cruising  and  high  speeds.  The  value  of  — became 
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less  positive  as  the  angle  of  attack  was  increased  in  all 
cases  where  it  was  of  appreciable  magnitude  but  in 
no  case  became  negative.  It  therefore  appears  that 
the  parasite  pitching  moment  will  be,  in  general,  such 
as  to  decrease  —  nmn  with  corresponding  effects  upon 
the  stability  characteristics. 

The  tail  length,  size,  and  aspect  ratio  are  very 
important  factors  in  the  stability  characteristics.  Tail 
length  appears  both  in  —mq  and  in  It  will  be 

remembered  that 


—  m^—Kv 


'2 ky2  S  da, 


m  Ic  dCL.n 
~*m°-  pSl  2kr\_da  {C"' 
da.t  dCZf  gt 


-0,)  + 


dCx  z 


dC\ 


da  da 


da  da(  S 


>3 


Increasing  l  has  a  very  beneficial  effect  upon  the  factor 
—  mqi  which  is  proportional  to  the  square  of  l.  In 
an  actual  case  increasing  l  will  also  increase  kYi  and 
perhaps  in  some  cases  rftt  so  that  the  net  increase  in 
— mq  will  only  approximate  that  indicated  by  consider¬ 
ing  only  the  square  of  the  tail  length.  Increasing  l 
makes  — larger  because  of  the  increase  in  moment 
arm  of  the  tail.  There  are  also  secondary  effects  because 


l  affects  -gjS  t iu  and  kY .  The  net  effect  of  small  in¬ 
creases  in  tail  length  can  only  be  predicted  from  con¬ 
sideration  of  the  particular  design.  There  is  a  small 
range  of  conditions,  corresponding  to  small  values  of 
— }imat  in  which  increasing  l  increase's  the  time  to  damp 
because  the  effect  upon  —  nma  is  of  more  importance 
than  the  increase  in  —  mq.  Large  increases  in  l  will 
always  increase  the  period  and  decrease  the  time  to 
damp. 


It  has  been  mentioned  that  l  affects  Tjt 


and 


da^ 

da 


The  effect  upon  rj,  is  small  and  may  he  neglected. 

da 

Increasing  tail  length  increases  The  angle  of 


attack  at  the  tail  may  be  expressed  as 


ai~a-\~it  —  € 

where  it  is  the  incidence  of  the  tail  referred  to  the  wing 
chord. 

The  downwash  may  be  expressed  as 


«”^(x+1)“0M(2/+1)~°‘23Ciz,  (reference  10) 
ti 


where  R  is  the  equivalent  monoplane  aspect  ratio  of 
the  wing  combination  and  x  and  ?/  are  the  distances  of 
the  tail  plane  behind  and  above  or  below  the  trailing 
edge  of  the  wing  in  chord  lengths.  The  value  e 
therefore  decreases  exponentially  with  tail  length  and 
hence 

dat  ^ _ dt 

da  da 


=  l-K^ 

da 


becomes  greater  with  increase  of  !. 


The  size  of  the  tail  appears  in  the  ratio  SJS  occurrin 
in  the  expressions  for  —  mq  or—  m«.  Increasing  Stj 
increases  both  —  mq  and  In  the  case  of  —n 

the  increase  is  proportional  to  St/S  but  in  the  cut 
of  the  increase  is  more  than  proportional  I 

the  increase  in  SJS  because  ~txtna  would  normally  1 
positive  were  SJS  equal  to  zero.  Increasing  tail  si? 
may  be  expected  to  have  a  net  effect  somewhat  les 
beneficial  than  increasing  /  unless  the  airplane  alread 
has  a  large  value  of  —  }±ma. 

The  aspect  ratio  of  the  tail  surfaces  affects  tl 


stability  through 


dCZt 

do i( 


Increases  in 


dCzt  .  , 

—  have  the s am 

da  t 


effect  as  increases  in  jStfS  and  need  no  further  disem 
sion.  It  will  be  noticed  that  if  SJS  is  increased  b 


dC 

increasing  the  chord,  — ~  will  bo  decreased  so  that  th 

dat 

net  effect  will  be  less  than  would  be  expected  on  th 
basis  of  tail  size  alone. 

The  tail  efficiency  factor  ijt  affects  the  stabilit 
characteristics  in  the  same  manner  as  the  ratio  SJt 
Tail  efficiency  depends  upon  the  interference  effects  o 
the  wings,  landing  gear,  and  fuselage  upon  the  ta: 
plane.  Warner  states  (reference  11)  that,  in  general 
tail  efficiency  has  been  found  to  range  from  0.6  to  0. 
with  extreme  values  as  low  as  0.54  and  as  high  as  l.U 
Computations  made  by  the  author  for  a  number  o 
military  airplanes  indicate  that  for  modern  design 
the  tail  efficiency  factor  will  be  of  the  order  of  0.75  t< 


0.8. 


The  series  of  charts  of  stability  characteristics  agains 

—  mq  and  —tima  has  been  extended  to  include  th 
range  of  high-lift  devices  known  at  the  present  tirm 
(figs.  41  to  54).  The  charts  bear  no  evidence  of  in 
creased  difficulty  with  longitudinal  stability  by  th< 
use  of  those  devices.  Devices  which  increase  th 
parasite  drag  as  well  as  the  maximum  lift  eoefficien 
are  more  desirable  than  those  which  merely  serve  (< 
maintain  the  flow  to  high  angles  of  attack.  It  wil 
be  noticed,  however,  that  in  order  to  use  the  chart 
with  variable-area  lifting  surfaces  the  true  area  of  tin 
extended  wing  must  be  used  in  order  to  come  withii 
the  range  of  parameters  given.  When  the  true  are; 
is  used,  SJS  and  m  will  have  been  decreased  and  th< 
quantity  c{Cg—Ca)  changed.  The  effective  value  o 

—  mq  is  thereby  decreased  as  compared  with  the  win<. 
with  lifting  surfaces  contracted;  whereas  the  value  o 

—  tima  may  be  either  decreased  or  increased,  depending 
upon  the  relative  values  of  the  quantities.  The  far1 
that  an  airplane  has  good  stability  characteristics 
with  conventional  surfaces  is  therefore  no  assurance 
that  trouble  will  not  be  experienced  when  high-lifi 
devices,  particularly  those  increasing  the  wing  area 
are  added.  In  connection  with  calculations  of  stability 
with  high-lift  devices  it  should  be  noted  that  in  cases 
of  such  devices  as  split  flaps,  variable-area  wings,  etc., 
which  greatly  increase  the  camber  of  the  wings,  par¬ 
ticular  attention  should  be  paid  to  reductions  in  tail 
efficiency  which  may  occur  in  some  instances 
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because  of  the  change  of  effective  location  of  the  tail 
surfaces  relative  to  the  wing  chord.  (See  also  refer¬ 
ence  12.) 

The  analysis  and  discussion  have  been  confined  to 
consideration  of  only  the  long-period  oscillations  or 
the  divergence  of  an  airplane  as  a  rigid  body.  The 
effects  of  such  physical  characteristics  as  wing  or 
tail-supporting  structure  elasticity  or  of  free  longitudi¬ 
nal  controls  are  not,  however,  unimportant.  No 
attempt  will  be  made  to  analyze  these  effects  mathe¬ 
matically  but  they  will  be  briefly  discussed  in  connec¬ 
tion  with  the  basic  relationships  brought  out  in  the 
foregoing  portions  of  the  report. 

The  effect  of  wing  elasticity  has  been  treated  mathe¬ 
matically  in  reference  13.  Although  these  effects  are 
normally  not  important  they  may  in  some  cases, 
particularly  in  view  of  the  high  speeds  being  obtained 
with  modern  aircraft,  be  sufficient  to  change  a  nor¬ 
mally  positive  value  of  —nma  to  a  negative  value  with 
possibly  serious  tendencies  to  divergence.  The  obvi¬ 
ous  means  of  guarding  against  such  possibilities  are 
to  make  the  wings  rigid  in  torsion  and  to  use  airfoil 
sections  having  positive  or  zero  values  of 

The  effect  of  elasticity  of  the  tail-supporting  struc¬ 
ture  is  primarily  of  importance  in  connection  with  the 
phenomena  of  buffeting  and  tail  flutter  and  needs  no 
consideration  here  other  than  to  call  attention  to  the 
possibility  of  resonance  between  the  natural  frequency 
of  oscillation  of  the  tail  structure  and  the  short-period 
oscillations  of  the  airplane,  which  are  normally  very 
heavily  damped.  The  results  of  such  resonance  would 
probably  be  structural  failure  and  would  not  receive 
practical  classification  as  instability. 

The  effect  of  freedom  of  the  longitudinal  controls 
has  been  treated  mathematically  in  reference  14.  The 
mathematical  analysis  introduces  several  new  factors 
into  an  already  complex  problem  and  it  seems  likely 
that  no  designer  would  attempt  to  predict  the  stability 
with  free  controls.  From  a  practical  standpoint  it  is 
obvious  that  if  the  control  surfaces  do  not  deflect 
appreciably  as  the  airplane  follows  the  sinusoidal  vari¬ 
ation  of  angle  of  attack,  etc.,  the  stability  will  be  the 
mimo  with  free  controls  as  with  fixed  controls. 

The  influences  that  cause  the  control  surfaces  to 
deflect  are  the  aerodynamic  hinge  moment  caused  by 
change  of  angle  of  attack  of  the  tail  and  air  speed ;  the 
inertia  moment  that  tends  to  keep  the  surfaces  from 
following  the  changes  of  attitude  of  the  airplane:  and 
the  mass  moment,  in  the  case  of  statically  unbalanced 
controls,  that  tends  to  deflect  the  elevator  under  the 
influence  of  gravity  and  accelerations  of  the  tail  plane 
is  accelerated  normal  to  the  longitudinal  axis.  The 
aerodynamic  hinge  moment,  which  is  normally  of  such 

sign  as  to  decrease  and  hence  both  -  m/zff  and  —  vmai 

may  be  kept  small  by  aerodynamic  balance  of  the 
elevator  surfaces.  The  inertia  moment  is  small  and 


may  be  neglected,  as  in  reference  14.  The  mass 
moment  may  be  kept  small  by  statically  balancing 
the  elevator. 

It  follows  from  the  foregoing  discussion  that  sta¬ 
bility  with  free  controls  can  be  insured  bv  providing 
ample  stability  with  fixed  controls  and  then  aero- 
dynamically  and  statically  balancing  the  elevator  sur¬ 
faces  or  making  the  control  irreversible.  It  is  generally 
not  desirable  to  attain  complete  aerodynamic  balance, 
but  the  difference  between  the  stability  with  free  con¬ 
trols  and  that  with  fixed  controls  can  be  made  small 
without  complete  aerodynamic  balance  or  without 
making  the  controls  completely  irreversible. 

CONCLUSIONS 

1.  It  is  believed  that  for  the  average  student  of  the 
subject  this  substitute  treatment  of  the  derivation  of 
the  relationships  governing  longitudinal  stability  will 
be  more  understandable  than  the  classical  treatment. 

2.  Increase  in  the  amplitude  of  longitudinal  oscilla¬ 
tions  with  time  is  possible  only  when  the  change  in 
angle  of  attack  lags  the  change  in  flight-path  angle. 
The  lag  of  the  angle  of  attack  with  respect  to  the 
flight-path  angle  depends  primarily  upon  the  static- 
stability  factor  ~pLfna. 

3.  The  nondime nsional  parameter  —  nma  not  only 
plays  the  mo3t  important  part  in  determining  the  lag 
of  the  angle  of  attack  with  respect  to  the  flight-path 
angle  but  also  is  the  only  source  of  unstabilizing 
angular  momentum  with  power  off.  Increasing  the 
value  of  —nma,  i.  e.,  increasing  the  static  stability, 
will  increase  the  tendency  to  instability  if  the  original 
value  of  — fjLma  is  small.  After  a  certain  value  of 
—  fjLma  is  reached,  further  increases  in  the  negative 
sense  will  reduce  the  tendency  to  instability. 

4.  The  effect  of  the  verticil  location  of  the  wing 
relative  to  the  center  of  gravity  should  not  be  neglected 
in  calculating  the  static-  stability. 

5.  Aerodynamically  efficient  airplanes  are  more 
likely  to  be  dynamically  unstable  than  less  efficient 
types. 

6.  Increasing  the  wing  loading  will  increase  the 
period  of  the  oscillations  but  will  not  necessarily  make 
the  airplane  either  more  or  less  stable. 

7.  High-lift  devices  that  depend  partly  upon  in¬ 
creases  in  wing  area  may  lead  to  instability  if  the  size 
and  disposition  of  the  tail  surfaces  are  determined  on 
the  basis  of  the  contracted  wing. 

8.  There  is  no  likelihood  of  dynamic  instability  at 
the  angle  of  attack  corresponding  to  maximum  or 
cruising  speeds.  The  damping  of  the  oscillations  may 
be  more  or  less  pronounced  at  low  speeds  than  at  high 
speeds.  There  is  likelihood  of  dynamic  instability  at 
low  speeds. 

9.  Stability  with  free  controls  can  be  brought  close 
to  that  with  fixed  controls  by  aerodynamically  and 
statically  balancing  the  elevators. 


Rotational  damping  factor ,  -  rn^  Rotational  damping  Factor, 
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Figures  15  to  20.— Dynamic  stability  charts. 
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Figures  21  to  26.— Dynamic  stability  charts. 
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Figures  27  to  32.— -Dynamic  stability  charts. 
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Figures  45  to  50. — Dynamic  stability  charts. 
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FIGUE*  63. 


Figures  51  TO  54.— Dynamic  stability  charts. 


FIGURE  64. 
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II.  ESTIMATION  OF  STABILITY 

METHOD 

The  charts  included  in  this  report  (figs.  15  to  54) 
provide  an  easy  method  of  estimating  the  dynamic 
stability  of  airplanes.  In  order  to  use  the  charts  it  is 
necessary  to  know  the  values  of  the  nondimensional 

parameters  CLi  CDt  —  aIK*  —  ^ma  ^or 

airplane  at  the  speed  under  consideration.  Knowing 
the  parameters,  turn  to  the  charts  and  select  those  for 

which  CLt  CDi  ^  and  most  neariY  represent  the 

case  under  consideration.  Spot  the  particular  point 
representing  —  and — y.vna  upon  the  charts.  Rough 
interpolations  between  the  values  of  f'  and  so  in¬ 
dicated  will,  in  general,  provide  a  sufficiently  accurate 
estimation  of  the  stability  characteristics  of  the  design. 
The  time  to  damp  to  one-half  amplitude,  in  seconds, 


T= 


-0.313 


VP 


r' 


and  the  period  in  seconds  is 

r  „  2.83V  s  L 

L  *' 

The  number  of  oscillations  for  the  amplitude  to  die  to 

T  V 

one-half  the  original  value  will  be  p=  —  g-p* 

DETERMINATION  OF  NONDIMENSIONAL 
PARAMETERS 


The  nondimensional  parameters  should  be  taken 
from  available  flight  or  wind-tunnel  data  if  possible 
If  such  data  cannot  be  obtained  they  may  be  calculated 
by  methods  outlined  in  the  following  paragraphs. 

Lift  coefficient  CL. — The  lift  coefficient  CL  follows 
directly  from  the  wing  loading  and  the  speed.  In  gen¬ 
eral,  it  is  desirable  to  choose  CL  as  0.20,  0.80,  1.40, 
1.90,  or  2.40  to  facilitate  use  of  the  charts.  For  air¬ 
planes  with  variable-area  wings  the  wing  area  in  use 
under  the  flight  condition  should  be  used  as  the  basis 
of  calculations. 

Drag  coefficient  CD.— The  drag  coefficient  CD  may 
be  estimated  from  the  relationship 


(\>  (\> 


CT, 2 

{kb)2 

r—che 

o 


(see  reference  (V) 


where  CDp  is  the  parasite-drag  coefficient, 


(reference  7  gives  k  in  convenient  form). 
e,  an  airplane  efficiency  factor  ranging  from  1 
for  very  clean  designs  to  0.7  for  very 
inefficient  designs  ('reference  6  V 


Slope  of  lift  curve  — The  slope  of  the  lift 

curve  may  be  taken  as  4  (all  values  of  angles  are 

in  radian  measure)  at  low  and  intermediate  angles  of 
attack  and  as  3  near  the  stall  if  the  values  of  the  other 
parameters  are  based  on  the  true  wing  area  and  the 
airplane  is  a  conventional  type.  It  should  be  noted 
that  in  cases  of  extremely  high  or  extremely  low  aspect 
dC 

ratios  for  which  differs  considerably  from  4  the 

wing  area  may  be  considered  such  as  to  make 

^=4  if  CL,  CD}  and  - m n,  and  -ma  are 

each  multiplied  by  the  ratio  of  4  to  the  actual  value  of 

for  reference  to  the  charts.  In  such  a  case  the 
da 

ratio  must  not  be  neglected  in  the  conversion  to  time 
in  seconds. 

Slope  of  the  drag  curve  — The  slope  of  the 

drag  curve  may  be  estimated  from  the  relationship 


da  ”  (kb)2 


The  value  of  here  used  may  be  taken  as  4  through¬ 
out  the  flight  range  as  ^  does  not  fall  off  with  ^ 


near  the  stall. 

Rotational  damping  factor  —■ mq. — This  parameter 
will,  in  general,  have  to  be  estimated  although  it  may 
be  determined  from  tests  of  an  oscillating  model,  or 
from  tests  on  a  whirling  arm  (reference  15).  For 
purposes  of  estimation  the  equation 

l2  St  dCLt 

~-mq  —  ixru  2kyl  s  dat 

where  K  is  taken  as  1.25  will  generally  be  sufficiently 
accurate.  The  tail  efficiency  t u  may  be  taken  as  0.75 
or  0.80  for  modern  designs.  The  distance  from  the 
center  of  gravity  to  the  elevator  quarter-chord  point 
l,  the  ratio  of  tail-plane  area  to  wing  area  St/S,  and 
the  radius  of  gyration  of  the  airplane  about  the  span 
axis  kY  are  dimensional  characteristics  of  the  airplane. 
The  radius  of  gyration  may  be  found  from  the  rela¬ 
tionship 


where 

B  is  the  moment  of  inertia  about  the  span  axis. 
g ,  the  acceleration  of  gravity. 

W ,  weight  of  the  airplane, 
or  estimated  from  the  relationship 

ky=^!CM2+W) 


28 


REPORT  NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


where 

CB  is  an  empirical  constant. 

/i,  over-all  length  of  the  airplane. 
h,  over-all  height  of  the  airplane. 

From  11  airplanes  for  which  CR  was  determined 
(reference  16) 

0.0325^  Cb^.  0.0394  with  the  average  value  0.0362. 
In  the  absence  of  test  data  the  slope  of  the  lift 
curve  for  the  tail  may  be  taken  as 
dCLl  ,5  5 

~~j — = — -r-  (reference  17) 


Static  stability  factor 
factor  is 


-The  static  stability 


It  dCm 


<S’pX32.2X/  2 kY2  da 


W 

6.5^c  dCm 


- j- j - ^ —  (standard  conditions) 

where  c  is  the  chord  length  upon  which  Cm  is  based. 
Various  methods  have  been  proposed  for  estimating 
dC 

■  g  •  The  author  prefers  the  relationship 

JCL(„  dCXtoz_  dCnv  l  St  dCLt  dat 

da  da''  9  a'  da  c  da  S  dat  da 

where 

C0  is  the  distance  of  the  center  of  gravity  of  the 
airplane  from  the  leading  edge  of  the  reference 
chord  in  chord  lengths. 

Caf  the  distance  of  the  aerodynamic  center  of  the 
wings  (reference  8)  from  the  leading  edge  of 
the  reference  chord  in  chord  lengths. 

dCXw 

~~l~'  the  slope  of  the  curve  of  longitudinal  force  co¬ 
efficient  for  the  wings  alone  against  angle  of 
attack  and  is 

d€xw  / '  dCDu\  /  dCL\  . 

~ir=\CL~~^r ) cos  sm  a 

dCr> 

where  and  CDw  are  values  for  the  wing 

alone. 


The  value  of  Ca  varies  from  0.23  to  0.25  refen 
to  the  mean  aerodynamic  chord  of  conventional  v. 

dC 

foils  (reference  9).  The  value  of  ^ p  may  be  e^ 

mated  if  the  designer  so  desires.  The  author  1 
found  the  value  to  vary  from  0  to  0.4  for  a  few  milita 
airplanes  at  low  angles  of  attack  with  the  avera 

value  about  0.2.  The  value  of  -  was  found,  i 

da 

the  cases  considered,  to  fall  off  to  approximately  z< 
at  high  angles  of  attack. 

The  rate  of  change  of  downwash  angle  at  the  t 
with  change  of  angle  of  attack  may  be  estimated  fn 
the  relationship 


£=m(x+1)' 


%+n-°- 


3<rfa  (re^erence  10) 


is  in  radian  measure 
da 

x}  the  distance  of  the  tail  plane  to  the  rear  of  tl 
trailing  edge  of  the  wing  in  chord  lengths. 
y,  the  distance  of  the  tail  plane  above  or  belo 
the  trailing  edge  of  the  wing  in  chord  length 

For  the  average  case  y— 0,  £=2.5, 

~~S~ 

EXAMPLES 


(7D=0.080  ^=0.39 

aa 

(Taken  from  full-scale  flight  data.) 

—  mq= 2.6  — /umtt=16.5 

(Estimated.) 

From  figure  23  f'=— 0.44  ^'=0.48 

From  figure  25  f'  =  —  0.035  i//  =  0.48 
From  interpolation  between  the  figures  on  the  bas; 
of  CD,  £'=  —  0.039  and  ^'  =  0.48.  From  interpolatio 


and  CDw  are  values  for  the  wing  between  the  figures  on  the  basis  of  ■ 


— »  the  distance  of  the  reference  chord  below  the 

c 

center  of  gravity  in  chord  lengths. 

dCmo 

the  slope  of  parasite  pitching  moment  due  to 

fuselage,  engine  nacelles,  landing  gear,  etc., 
against  angle  of  attack. 

the  rate,  of  change  of  the  angle  of  attack  at  the 

tail  with  change  of  angle  of  attack  of  the  air¬ 
plane  and  is 

j^=  1  — ^  where  e  is  the  angle  of  downwash  at  the  tail. 


and  ^  =  0.48.  From  full-scale  unpublished  flight  test 
of  this  airplane  £'=  —  0.041  and  ^'  =  0.72. 


CD= 0.076  ^=0.42 

da 

—mq=  1.5  —^ma— 9 

(Taken  from  full-scale  flight  data.) 

From  figure  23  f'  =  -0.034  ^'=0.49 

From  figure  24  {*'=— 0.030  i//— 0.50  - 

From  figure  25  £'  -  —  0.025  ^'—0.48 

From  interpolation  between  figures  23  and  24  on  th< 
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basis  of  CD,  £'  =  -0.031  and  *'=0.50.  From  inter¬ 
polation  between  figures  23  and  24  on  the  basis  of 

<i£pi  —  0.033  and  *'=0.49.  From  interpolation 
da  ^ 

between  figures  24  and  25  on  the  basis  of 

£'  =  —0.027  and  *'  =  0.49. 

From  full-scale  flight  tests  of  this  airplane  (refer¬ 
ence  2)  £'  =  —  0.038  and  *'  =  0.55. 

AIRPLANE  C 

Ci=0.80  ^=4.00 

<70=0.094  ^=0.51 

~m?=2.3  —nma—4. 8 

(Estimated  data.) 

From  figure  23,  £'  =  —0.034  and  *'=0.4. 


No  interpolation  is  necessary.  From  extrapolation  of 
full-scale  unpublished  flight  data  £'  =  — 0.030  and  *'  = 
0.72. 

AIRPLANE  D 

<7t=0.80  ^=4.00 

Co- 0.12  ^>=0.54 

—  mQ~  1.76  — ^ma  =  3.40 

(Estimated  data.) 

From  figure  26,  £'  =  —0.03  to  —0.04  and  *'  =  0.39. 
From  full-scale  unpublished  data  £'  =  —0.035  and  *'  = 
0.53. 


Langley  Memorial  Aeronautical  Laboratory, 
National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  December  13,  1934 


APPENDIX 


A, 

B, 

C, 

D, 

E, 


RECAPITULATION  OF  SYMBOLS 


^coefficients  of  stability  quartic. 


a* | coefficients  of  stability  quadratic. 

c=K3C,°+© 


c^+c*) 


Hc«2 

e,  base  of  natural  logarithms. 

e,  airplane  efficiency  factor  such  that  Cor- 


Ql 

Akbf 


Cr ,  coefficient  of  resultant  force. 

CZ)  coefficient  of  force  normal  to  the  wing  chord, 
positive  downward. 

Cxt  coefficient  of  force  parallel  to  the  wing  chord, 
positive  forward. 

Ca,  coefficient  of  aerodynamic  center  (ratio  of  dis¬ 
tance  of  aerodynamic  center  from  leading  edge 
to  chord  length). 

Cgy  coefficient  of  center  of  gravity  (ratio  of  distance 
of  center  of  gravity  from  leading  edge  to  chord 
length). 

Cmp,  coefficient  of  pitching  moment  other  than  that 
from  wings  and  horizontal  tail  surfaces. 

Cm0,  coefficient  of  pitching  moment  of  wing  at  zero 
lift. 

Cmc  g ,  coefficient  of  pitching  moment  with  respect  to 
center  of  gravity. 

C  —  frr2  coefficient  of  radius  of  gyration  about 

B~~ {li2-\-h2y  lateral  axis. 

llt  over-all  length  of  airplane. 

ft,  over-all  height  of  airplane. 

/,  distance  from  center  of  gravity  to  quarter-chord 
point  of  horizontal  tail  surfaces. 

i  =  -yj —  1 

K,  empirical  factor  by  which  computed  value  of 
rotational  damping  factor  of  tail  is  multiplied 
to  give  total  rotational  damping  factor. 

ft,  empirical  factor  for  conversion  of  biplane  span 
to  equivalent  monoplane  span. 
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lc  dCmc  0 

(HT 

l2  St  dCL( 

2 ft?  S  ~d^ 7* 

0,  subscript  denoting  that  the  value  of  the  symbol 
refers  to  zero  time. 


m 

2 Tc? 


—  ma=K  i 


->  static-stability  factor, 
rotational  damping  factor. 


P',  period  of  an  oscillation,  nondimensional  units. 

P,  period  of  an  oscillation,  seconds. 

T' ,  time  for  an  oscillation  to  decrease  to  one-half 
amplitude,  nondimensional  units, 

T,  time  for  an  oscillation  to  decrease  to  one-half 
amplitude,  seconds. 

t,  subscript  denoting  that  the  symbol  refers  to  the 
horizontal  tail  surfaces. 

&V',  ratio  of  A V  to  V. 

w ,  subscript  denoting  that  the  symbol  refers  to  the 

wing  alone. 

x,  distance  of  horizontal  tail  surfaces  to  the  rear 

of  the  trailing  edge  of  the  wing,  chord  lengths. 

y ,  distance  of  horizontal  tail  surfaces  above  or 

below  the  wing  chord  (extended),  chord 
lengths. 

2,  distance  of  wing  chord  below'  center  of  gravity. 
$i,  angle  of  lag  of  change  in  angle  of  attack  wuth 
respect  to  change  in  angle  of  flight  path. 

82,  angle  of  lag  of  change  in  velocity  along  flight 
path  with  respect  to  change  in  angle  of  flight 
path. 

53,  angle  of  lag  of  change  in  angle  of  pitch  with 
respect  to  change  in  angle  of  flight  path. 

8{,  angle  of  lag  of  change  of  angle  of  attack  with 
respect  to  change  in  angle  of  pitch. 
rjt,  tail  efficiency. 

t'  =  —■ — damping  coefficient. 


=  period  coefficient. 

x'=r'±tV'. 

n  =  relative  density. 

pot 

r—  — time  conversion  factor, 
po  V 
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4.  PROPELLER  SYMBOLS 


D,  Diameter 

p,  Geometrio  pitch 

p!D,  Pitch  ratio 

V',  Inflow  velocity 

V„  Slipstream  velocity 

T 

T,  Thrust,  absolute  coefficient  6V 

Q,  Torque,  absolute  coefficient  CQ = 


Power,  absolute  coefficient  <?*>=•  — 

Speed-power  coefficient  . 

Efficiency 

Revolutions  per  second,  r.p.s. 
Effective  helix  angle  =  tan-1 


*  ~  “  0  ?  f  K 


5.  NUMERICAL  RELATIONS 

i  TK  —  n  i.-< 


1  metric  horsepower  = 1.0132  hp. 
1  m.p.h.= 0.4470  m.p.s. 

1  m.p.s.  =2.2369  m.p.h 


1  kg  =  2.2046  lb. 

1  mi.  =  1,609.35  m  =  5,280  ft. 
lm  =  3.2808  ft. 


